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Abstract 

The fate of symmetries at high temperature determines the dynamics 
of the very early universe. It is conceivable that temperature effects favor 
symmetry breaking instead of restoration. Concerning global symmetries, 
the non-linear sigma model is analyzed in detail. For spontaneously broken 
gauge symmetries, we propose the gauge boson magnetic mass as a "flag" 
for symmetry (non)-restoration. We consider several cases: the standard 
model with one and two Higgs doublets in the perturbative regime, and 
the case of a strongly interacting Higgs sector. The latter is done in a 
model independent way with the tools provided by chiral Lagrangians. Our 
results clearly point towards restoration, a pattern consistent with recent 
lattice computations for global symmetries. In addition, we explicitly verify 
BRST invariance for gauge theories at finite temperature. 



LPTHE Orsay-98/02 
FTU AM/97/20 
FTUV/97-73 
IFIC/97-105 
hep-ph/ 980 1244 



1 Laboratoire associe au Centre National de la Recherche Scientifique - URA D00063 
e-mail: gavela@delta.ft.uam.es, nuria@goya.ific.uv.es, pene@qcd.th.u-psud.fr, 
vargas@delta.ft.uam.es. 



2 



In which sense does one say that an internal symmetry is restored or broken 
due to temperature effects? What is the relevant order parameter? And whenever 
more than one such parameter can be defined, for which physical consequences are 
their differences relevant? These are the type of questions to face when discussing 
symmetry (non-)restorationF]. 

The vacuum structure of a system remains unchanged when it is heated. 
In this sense the degree of symmetry of a system is not modified. "Symmetry 
restoration" due to temperature effects is thus a misleading denomination for a 
very simple effect: the spontaneous breaking of a global or gauge symmetry can 
be masked for all physical purposes when thermal agitation is present. This suits 
intuition, as a thermal excitation gives in general a positive energy contribution, 
allowing particles to "climb" barriers between separate minima and finally hiding 
those barriers for high enough temperatures. Thermal field theory computes 
these effects and usually synthesizes them in the form of a so-called effective 
potential whose minimum sits at zero values of the fields. Ferromagnets provide 
well-known experimental examples of a similar behavior when heated above some 
critical temperature. 

The suggestion that spontaneously broken field theories are restored at high 
temperature was first made by Kirzhnits and Linde They gave qualitative 
arguments to support this idea in the case of global symmetries. In the same 
direction pointed the results of Dolan and Jackiw and Weinberg [[| for gauge 
theories (although in this case the choice of the scalar field vacuum expectation 
value as an order parameter is a delicate one). 

Weinberg noticed as well an opposite possibility: global symmetry non-res- 
toration at high temperatures for scalar potentials with more than one Higgs 
multiplet. With just one Higgs the scenario is ruled out due to the constraints 
imposed on the scalar self-coupling by the boundedness of the potential, while 
models with two (or more) multiplets can easily accommodate it. The same be- 
havior was found in the Schwinger model and in a dynamical model of symmetry 
violation in four dimensions H. 

An analogous situation has been experimentally observed in nature for the 
ferroelectric material known as Rochelle salt, which shifts from a disordered phase 
to a more ordered one when heated, as measured by the spontaneous polarization 
parameter. In the case of the Rochelle salt the symmetry is restored again for 
high enough temperatures, though. Common sense suggests that this should be as 
well the case in field theory, with thermal excitations dominating the free energy 
unless some finite parameter, such as finite volume, causal domain size, etc., plays 
a role. Without entering to discuss it, it is clear that even a temporal intermediate 
period, in which thermal effects enhance the effective symmetry breaking instead 

2 A related question is the so-called inverse symmetry breaking, describing systems for which 
the symmetry is exact at zero temperature and broken when heated; all through the paper we 
will take the liberty of dubbing symmetry non-restoration both scenarios, unless the contrary 
is explicitly stated. 
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of restoring it, could have far reaching cosmological consequences. 

It is worth to remark, though, that Weinberg results on symmetry non- 
restoration are based on the one-loop approximation to the finite temperature 
effective potential, which is known to be unreliable for the discussion of many 
aspects of phase transitions. Different techniques, including non-perturbative 
ones, are being actively applied to improve the one-loop approximation, mainly 
for the study of global symmetries. The results are very interesting and quite 
often contradictory: some studies confirm that symmetry non-restoration exists, 
although with a sizable reduction of the parameter space where it occurs 
while other analysis conclude that symmetry is always restored at high temper- 
ature when non-perturbative effects are taken into account It has been 
shown that in a finite lattice no order is possible at sufficiently high temperature 
nj. Although the relevance of this result for the continuum limit is unclear, a 



Monte Carlo simulation in 2+1 dimensions seems to support this conclusion [ITT 



Symmetry non-restoration is indeed being increasingly reconsidered as a can- 
didate way out of many cosmological problems arising in spontaneously bro- 
ken theories. Examples are the domain wall and axion problems [[12| and the 
monopole problem in Grand Unified Theories [[HI . 

As recalled in Section |3|, in the minimal standard SU(2) <g> £7(1) model the 
symmetry is necessarily restored, given the simplicity of its Higgs sector. At 
present, there are two main avenues to explore physics beyond the Standard 
Model: theories in which the Higgs particle is a fundamental one, supersymmetry 
being its most representative example, and those for which it is not, currently 
dubbed as strongly interacting Higgs scenarios. 

Supersymmetry is broken "de facto" at high temperatures, due to the differ- 
ence in the boson and fermion populations, as dictated by Bose-Einstein versus 
Fermi-Dirac statistics. The debatable and interesting question is whether the 
internal symmetries present in supersymmetric theories, and whose fate is funda- 
mental for the existence of topological defects, are restored. It has been proven 
that such is the case for renormalizable supersymmetric theories WM. For the 



latter, a recent analysis for systems involving non-vanishing background charges 



shows that symmetry non-restoration could be possible [I3|. The consideration 
of non-renormalizable terms in the Lagrangians has led as well to a polemics: 
their mere addition does not lead to symmetry non-restoration . 

Here we rather follow the path leading to a non-elementary Higgs scenario. 
In so doing, we first reanalyze the global SU(Nf)n ® SU(Nf)i non-linear sigma 
model, relevant in supergravity and many other scenarios, in sect. [| Sect. |3] is 
devoted to the analysis of gauge symmetries; after discussing BRST invariance 
at finite temperature, we study the behavior of the 577(2) ® C7(l) symmetry in 
several scenarios. In subsect. [3.3| we analyze both the minimal standard model 
and the standard model with two Higgs doublets within the perturbative regime, 
while in subsect. we consider a strongly interacting Higgs sector in a model 
independent way, using the techniques of chiral Lagrangians, and we discuss the 



4 



differences with the results of the previous section. These different chapters are 
preceded by some comments on order parameters, sect.[j], and followed by our 
conclusions. 



1 The order parameter 

The interesting order parameter to consider in a phase transition depends first 
of all on the question one wants to study. 

An illustrative example is provided by spin systems in solid state physics. 
Both in ferromagnets and antiferromagnets, the ground state breaks rotational 
symmetry: the spins align for the former and display an antiparallel alignment 
for the latter. The traditional order parameter is the average spontaneous mag- 
netization (m) ^ which plays a crucial role in the description of the response of 
the system to an external magnetic field: it turns to be important for ferromag- 
nets, while marginal or even vanishing for antiferromagnets to the extent that the 
ground state approaches the Neel-type magnetic order. Hence, the spontaneous 
magnetization is an example of order parameter whose non-zero value is not nec- 
essary for the spontaneous breakdown of the symmetry. Ferrimagnets are yet 
another scenario: antialignment is present alike to the case of antiferromagnets 
although (m) ^ 0, as the weight allocated to the two possible spin projections 
differs. 

Analogous questions arise in particle physics: different so-called order param- 
eters can be correlated to different physical effects. The appropriate parameter 
depends on the aspect of the history of the universe under study, and not all of 
them necessarily "bip" simultaneously. 

Already at zero temperature, the relationship among different possible order 
parameters is not always straightforward. Recall massless QCD at low energies, 
with pion interactions appropriately described by chiral Lagrangians. The pion 
decay constant, F n , and the condensate, (^VP), are not necessarily equivalent 
order parameters. Although unnatural, = is not theoretically forbidden 

while a non-null v.e.v. of some higher dimension operator accompanies F n as a 



"flag" for dynamical symmetry breaking [JT7 . 

In a general way it is clear that when the Lagrangian, at zero temperature, is 
just a one parameter theory, all putative order parameters should be equivalent. 
Such is the case with most Lagrangians respecting global symmetries, where the 
value of the field at the minimum of the effective potential is a well-defined order 
parameter, commonly used, and any other one is simply related to it. 

On the contrary, for spontaneously broken gauge theories the issue is much 
more subtle. To begin with, the v.e.v. of any non-gauge invariant operator is 
necessarily zero 113. Only gauge invariant operators, such as |0| 2 , may have a 



non- vanishing v.e.v., signaling the Higgs mechanism. Once a gauge-fixing proce- 
dure has been performed, a gauge-non invariant minimum of the effective poten- 
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tial may appear, which may be useful whenever its physical meaning is properly 
extracted in due respect of general Ward identities. The same applies to the 
v.e.v. of higher dimension operators. 

The effective potential itself is gauge dependent. However, the values of the 
effective potential at its local minima or maxima are gauge independent. Hence 
if there is a minimum of V((f>) with a value lower than V(0) in one gauge, then 
there will be such minimum in any gauge (although its position will generally be 
different) and the symmetries will definitely be broken. 

In practice, gauge- dependent correlation functions are often used in the study 
of phase transitions: the physical conclusions are expected to be rather close 
to those derived with gauge-invariant ones if the fluctuations of the scalar fields 
are small compared to their vacuum expectation values. This can be safe in 
the broken phase of the theory, while quite misleading in the symmetric phase, 
as recently discussed in ref. [[TJJ, where a detailed description of the zoo of 
correlation functions can be found as well. 

It is worth to briefly specify the "flags" for symmetry breaking discussed in 
the present paper: 

- for the global symmetries of the non-linear a model and its extensions in 
terms of chiral Lagrangians, we discuss both the pion decay constant F n and 
the vacuum expectation value of the condensate. The latter is defined from an 
effective potential and can be interpreted as the remnant of the disappeared sigma 
field. Both parameters are essentially equivalent since, at zero temperature, we 
are dealing with a one parameter theory. 

- for the gauge symmetry case, specifically the standard electroweak model 
and its extensions, we concentrate instead on particle masses. In the perturbative 
regime, both the negative Higgs "mass" and the magnetic mass for the gauge 
bosons are discussed. When the Higgs particle disappears from the spectrum 
and we enter the non-perturbative regime of the gauged non-linear sigma model, 
our order parameter will be the gauge boson magnetic mass. 

The magnetic mass squared is defined as the temperature dependent contribu- 
tion to the transverse part of the gauge boson self-energy, Hr(0, k), for vanishing 
three-momentum k. At the order we work it is gauge invariant. Notice that 
Weinberg || advocates the use of gauge invariant operators carrying moderate 
momenta and zero energy as order parameters. 

The analogous electric mass, whose square is defined by the longitudinal com- 
ponent of the gauge boson self-energy, n^(0, k) with k — > 0, is not a suitable 
parameter. Indeed, it tends to increase at high temperature even when the sym- 
metry is restored. The intuitive explanation is electric screening: some particles 
in the theory carry an electric charge. Already at one-loop order, thermal fluctu- 
ations pull charged pairs out of the vacuum to screen external charges. However, 
there are no fundamental particles in any gauge theory which carry a magnetic 
charge. Magnetic screening can presumably then only arise from non-perturbative 
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fluctuations which carry magnetic charge. 

A perturbative computation of the magnetic mass shows that it is exactly 
equal to zero at one loop in an unbroken gauge theory. For unbroken non-Abelian 
gauge theories, such as QCD, higher orders in perturbation theory suffer from 
infrared divergences, and a magnetic mass of order g 2 T is expected to be gen- 
erated non-perturbatively. In spontaneously broken gauge theories, such as the 
standard electroweak model and its extensions, no such divergences are present. 
Thus, we propose to use the magnetic mass as a "flag", expecting that even in 
perturbation theory it will show a tendency to vanish at high enough tempera- 
tures whenever symmetry restoration occurs. Of course, it will be a valid "flag" 
only when exploring the broken phase of the theory, for the reasons given above. 

Another pertinent point to recall is that temperature corrections break Lorentz 
invariance as the plasma sets a preferred reference frame. Assume for instance 
a zero temperature Lagrangian based on an internal symmetry. Certainly the 
mixed states describing the new "effective vacuum" may greatly differ from the 
real vacuum structure. What about the finite temperature effective Lagrangian 
itself? Up to which point its functional form may differ from the initial one? 
Non-zero temperature is tantamount to treat space and time differently: internal 
symmetries at the Lagrangian level, such as chiral symmetry, cannot be explicitly 
broken due to it. What is to be a priori expected is a splitting of any operator 
into its temporal and spatial components, with differing coefficients. For instance 
F v in the chiral non-linear sigma model will generate two different coupling con- 
stants at finite temperature, a temporal one, F*(T), and a spatial one, F£(T) pl|]. 
We leave the corresponding considerations for gauge theories for the beginning 
of sect. [| 

A necessary condition for symmetry restoration is that all possible order pa- 
rameters or "flags" for symmetry restoration do signal it. 

2 Global non-linear sigma model: the ef- 
fective Lagrangian 

The restoration of spontaneously broken global symmetries is discussed in this 
section within an effective Lagrangian approach. We consider the SU(Nf)n ® 
SU(Nf)i, non-linear sigma model, which may be defined by the Lagrangian 

C=~d lt Tr a & i TT a + ~d ll <T& t <T, (1) 

with the constraint 

Fl = a 2 + n 2 . (2) 

By convention, we take the scalar condensate in the direction of the a component; 
that is, at tree level 

<^ 2 > = Fl . (3) 
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Since the global symmetry is broken down to an SU(Nf) symmetry, there are a 
total of Nj — 1 Goldstone bosons, which are identified with the pion fields, jf. 
The constraint @ determines the a field in terms of the pion fields, so that in 
the non-linear Lagrangian only the latter appears. A non-linear redefinition of 
the fields is possible without changing the physical content of the theory, leading 
to different parametrizations. The so-called exponential representation can be 
described by the Lagrangian 

= I Fl Tr^Ud^), (4) 

where U is a SU(Nf) unitary matrix field 

U = exp (z ^) , (5) 

with T a the generators of SU(Nf), normalized as Tr (T a T b ) = 2<5 afe and [T a ,T b ] = 
2ifab c T c , being f a b c the structure constants of SU(Nf). 

As it is well known, all realizations of the non-linear chiral Lagrangian, such as 
the exponential one (f|), square root, Weinberg, etc. with Nf = 2(3), are low 
energy effective theories for QCD with 2(3) massless quarks, expressed in terms 
of Goldstone bosons and systematically expanded in powers of the Goldstone 
bosons momenta. As a consequence of the chiral symmetry, these models possess 
the remarkable property of universality: once the coupling constants have been 
adjusted (F n is the only one at lowest order) all physical predictions are the 
same. Therefore, the chiral Lagrangian not only parametrizes the dynamics of 
the Goldstone bosons that emerge in QCD but also of any other theory, such as 
the Higgs model, that follows the same symmetry breaking pattern. 

We have analyzed two order parameters: the pion decay constant F n , and the 
vacuum expectation value of the a field, (a). Notice that, while their behavior 
should be essentially equivalent, their precise variation rate with temperature 
may differ somewhat. Indeed, the constraint (fj) as a thermal average implies 
(a 2 ) = F 2 , while in general (a 2 ) ^ (o") 2 - Our treatment differs from previous 
ones in that we have considered them as Lagrangian parameters, whose variation 
with temperature is read from the one-loop effective Lagrangian we derive. 

Although we will just discuss below the calculation in the exponential rep- 
resentation, we have explicitly checked that the results of measurable quantities 
are the same in other parametrizations used in the literature, namely the square 
root and Weinberg representations. Of course, for quantities without a physical 
meaning, the temperature corrections can be representation dependent. 

We drop all temperature independent ultraviolet divergent quantities from our 
expressions, recalling that when a theory is renormalized at zero temperature no 
more infinities of that type appear at finite temperature. 
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Figure 1: One loop self-energy diagram for the pions. 



2.1 Temperature corrections to F n 



Temperature corrections to F n are obtained from an effective Lagrangian ap- 



proach. The chiral Lagrangian 
certain order, 



£(2) 



1 



6F2 



can be expanded in powers of (■§-) up to a 



((7r9 M 7r) (TT^T?) - (7f7r)9 M 7T^7f) + . . . (6) 



We have computed the one-loop temperature corrections to this Lagrangian 
to leading order T 2 /F 2 , and proved that they lead to an effective Lagrangian 
with the same structure as the tree-level one, albeit with two F^'s: a temporal 
one, F*, and a spatial one, F*. It encloses the full temperature effects in the 
renormalized (temperature dependent) parameters. 

Due to the derivative character of the interactions, a contribution to the 
kinetic energy term, at leading order T 2 /F 2 , is obtained when computing the one 
particle irreducible (1PI) two point Green function at one loop (Fig. 1). This 
term is absorbed by pion field renormalization. In the exponential representation 
used here, we find 

(N — 1)T 2 ' 



IT 



36F2 



(7) 



Diagrams in Fig. 2 contribute to the 1PI four point function at one loop, 
leading to different thermal corrections for the spatial and the temporal coupling 
constants in which F v splits at finite temperature, as mentioned above. To this 
order, it results 



Fl{T) 



Fjt 



(N - 1)T 2 

(N + 1)T 2 
24F2 



(8) 
(9) 



where N = Nj — 1 represents the number of pions. Both temperature dependent 
renormalized parameters F£(T) and F*(T) show a clear tendency to vanish at 
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Figure 2: One loop diagrams contributing to the 1PI four point Green function 
for the pions. 



high enough temperatures, pointing towards chiral symmetry restoration. We 
have as well explicitly checked that F£ so derived is representation independent. 

Thermal corrections to the pion decay constant have been computed in the 
literature following different approaches p2| , p3| , [f2l"H , |25|| , |26| . Our result for 
the effective spatial coupling, F*(T), is in agreement with those calculations of 
F n (T). In most of them, F n (T) is obtained from its usual definition (slightly 
modified at finite temperature through the two point function of the axial 
vector current, and there is no splitting between temporal and spatial couplings 
at one loop; it appears at two loops |[27|| . Notice that since we consider F n just 
as a parameter in the Lagrangian, it does not necessarily coincide with the pion 
decay constant as usually defined. 

To avoid technical complications, we have computed the pion field and F n 
renormalization from the lowest order terms in the expansion of the Lagrangian 

£(2) 

in powers of the pion fields; chiral symmetry ensures that all higher terms 
in the field expansion are consistently renormalized once 7? and F n have been 
renormalized from these lowest order terms. 



2.2 Temperature Corrections to the Condensate at One 
Loop 

Temperature corrections to (a) are computed through the addition of a small 
chirality breaking term, which makes the Lagrangian slightly asymmetric [], that 

3 We recall that the QCD scalar density if? if? whose vacuum expectation value represents the 
familiar QCD condensate, is equivalent to a = -g-Tr(U + U + ) since and Tr(U + U + ) can 
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is, we consider 

C = C^ + C B , (10) 

with 

C B = cv = ^Tr(U + rf). (11) 

2 

Expanding the last term in powers of j% in the exponential representation it is 
found 

C B = cF^ —nn + — (tttt) 2 + . . . (12) 

Following the effective Lagrangian approach, we compute the one-loop order T 2 
corrections to L through the 1PI zero, two and four point Green functions. The 
kind of diagrams involved in the calculation are vacuum energy ones for the zero 
point 1PI Green function and the same as in the previous section (see Figs. 1 and 
2), although with modified couplings, for the two and four point 1PI functions. 
Now, besides F n and the pion field, also the parameter c is renormalized. Since 
it only appears in the product cF n , there is an ambiguity, depending on which 
(spatial or temporal) F n (T) we consider, leading to 



c -m = c (i-^j, us) 

<W = °(i + 5D|)- < 14 > 

Again, the one-loop effective Lagrangian written in terms of the temperature- 
dependent parameters has the same structure as the tree-level one, space-time 
splitted, though. 

Notice that (minus) the first term in the expansion of C B , —cF n , can be 
interpreted as the vacuum energy density of the system. That is, the free energy 
of a system of free bosons, given by 

c s (T)F:(T) = c t (T)F t n (T). (15) 

Since the operator a can be obtained by deriving the bare Lagrangian with 
respect to the parameter c (see eq. (p])), we can also interpret the result as a 
thermal correction to the scalar condensate Q. Taking the derivative with respect 

be shown to transform in the same way under the chiral group. Lb is thus equivalent to a 
quark mass term. 

4 Recall that the thermal average of the operator a is defined as 

Wt = Tr(e-eB) (16) 

where Tre- pH = J[dU]e' " is the partition function and f3 = l/T. Thus, one can compute 
(ct)t as the derivative of the partition function with respect to the parameter c, at c = 0. 
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to (bare) c of the one loop effective Lagrangian, the temperature corrections to 
the condensate are found. Explicit chiral symmetry is recovered by fixing c = 
at the end of the computation. The final result is satisfactorily the same whether 
either the spatial set F£, c s or the temporal set F*, c* is used, leading to: 

Mt= W(l-iv^), (17) 



in agreement with [f24j], [p2]] , As can be seen in flT7j), the temperature cor- 
rection to the condensate also points towards chiral symmetry restoration. 

Notice that Lb is just the well known classical potential up to a minus sign. 
However, it was not possible to use the standard method for computing effective 
potentials [0 due to the presence of derivative couplings. Using a generalization 



of this method [24] the same result is recovered, as already mentioned. 



3 Gauge symmetry: SU{2) ® U{1) 

In this section we study theories with gauge group SU(2) <g> U{1). We consider 
the cases with one light Higgs doublet, two light Higgs doublets, and the generic 
one where the Higgs sector becomes strongly interacting, the latter done in a 
model independent way. Before entering into such details, we dwell again into 
the delicate issue of the "flag" for symmetry (non-)restoration for gauge theories, 
and in the fate of gauge invariance itself when a system is heated. 



3.1 The magnetic mass 

In a gauge theory, the pseudo-Goldstone bosons of the Lagrangian are unphysical 
fields, unlike the gauge bosons. 

As stated in sect. [1], we choose the gauge boson magnetic mass as our "flag" 
or indicator for symmetry (non-)restoration. 

At non-zero temperature, the self-energy tensor of the gauge boson may de- 
pend on the four- velocity of the plasma u^. Consequently, the gauge boson 
self-energy can be expressed as a linear combination of four possible tensors: g^ u , 
kuK, u^Uy and k^u v + k v u^. Some linear combinations of these tensors are usually 



chosen as the standard basis set [23], denoted A^, B^, C^ v and and defined 



in Appendix A. In this basis, the one loop gauge boson self-energy is written as 

W u = U T A^ + UlB^ + U D D^, (18) 

where the subscripts T and L denote transverse and longitudinal with respect to 
the spatial component k of the wave vector. 

This is analogous to the extraction of (4'5')t in QCD, by first adding an explicitly chiral 
symmetry breaking term m\E"J' to the bare Lagrangian, computing the temperature corrections, 
and deriving then with respect to m J22| . 
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The magnetic mass is denned as ITt(0, k), with vanishing k. At one- loop and 
leading order, O(gT), it will be shown to be gauge invariant both for the standard 
model and for its extensions considered below. The explicit computations will be 
focused in the W gauge boson mass. 

3.2 Checking gauge invariance: BRST identities 

Up to our knowledge, the Slavnov- Taylor identities at finite temperature have 
never been explicitly verified in the literature for the electroweak theory. We 
explicitly perform such a task in the present work, for the two-point functions of 
the theory. 

Indeed, one expects gauge invariance to be preserved at non-zero tempera- 
ture. A simple reasoning can be developed in the imaginary time formalism, 
where finite temperature just amounts to compactifying the time direction, that 
is, to perform a global "distortion" of the system. Gauge transformations are 
local ones by definition, and thus they should not be affected by global topolog- 
ical conditions. Once the gauge fixing procedure has been implemented, BRST 
invariance remains, and the corresponding Slavnov- Taylor identities are to be 
proven. 

One should realize that the proof is much more juicy than at zero temperature: 
there, quadratic divergences are disposed of by counterterms from the start, and 
the Slavnov- Taylor identities for such non-physical quadratically divergent terms 
are not even considered. At finite temperature those quadratic divergences are 
the source of the T 2 dependence. It is then mandatory, and new, to check the 
BRST identities on them. 

Both in the linear and the non-linear realizations of the SU(2) eg) U(l) gauge 
symmetry, the Ward identities relating the two point Green functions at one loop 
are given by 

k 2 (Tl% + 2M w Tl W7r± ) -M^IF* = 0, 
A; 2 (n| - 2iM z H Z7T °) - M 2 Z W° = 0, 

k 2 lP D = 0, 
k 2 (TP D z -iM z W*°) = 0, 
i!i z D - iM z ^H Z7T ° - IT = 0, 
£11^ + M w $L Wir± - II C± = 0, (19) 

where £ is the gauge fixing parameter in R^ gauges. IT^ , II f,, IV D and If^f are 
the form factors introduced in eq. ([18]) for the W — W , Z — Z, 7 — 7 and Z — 7 
self-energies, respectively. Yl Ww is defined from the two point Green function 
with external legs — as 

njf 7 ^ = k^U W7r± , (20) 
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and analogously for U Zw and Il 77r , while IF and Il c represent the charged and 
neutral Faddeev-Popov ghosts self-energies. 

3.3 Perturbative Higgs sector 
3.3.1 Minimal Standard Model 

Consider the minimal electroweak standard model, that is, with just one light 
Higgs doublet. Here, all the couplings of the theory are in the perturbative 
range, and we can rely on the one-loop approximation to the effective potential 
at non-zero temperature. 

It is well known that with just one Higgs doublet the gauge symmetry is always 
restored at high temperature. Indeed, given the simplicity of the potential, 

V(<j>) T=0 = -/i 2 (0V) + A(0W, (21) 

the condition that it has to be bounded from below forces the sign of A to be 
positive. The one-loop thermal corrections to the above potential can be readily 
computed in gauges by the usual methods @. In the high temperature limit 
(T ^> rrii, with the masses of all standard model particles) the leading order 
T 2 corrections are gauge invariant and read: 

5f = -g(6A + \g 2 + \g' 2 + 3h 2 + 3h 2 b + h 2 T ), (22) 
leading to the temperature dependent Higgs v.e.v. 

(23) 



T 2 

v(T) 2 = v 2 - — 



\ 3g 2 g' 2 h 2 h 2 h 2 
1 + — + — + — + — + — 
8A 8A 2A 2A 6A 



where v 2 = fi 2 /X denotes the Higgs v.e.v. at tree level and h t , hb and h T are the 
Yukawa coupling constants of the quarks t and b, and the r lepton, respectively. 

As expected, exactly the same behavior is seen from the Z and W gauge 
boson magnetic mass. The set of diagrams contributing to the Z self-energy at 
one loop are shown in Fig. |3]. The coupling constants are not renormalized at 
one loop (at order T 2 ), which allows to write 

M 2 w , mag = 9 ^v{Tf, (24) 



Z,mag ~ ^ V (T) ? 



T 2 



, 3g 2 g' 2 h 2 h 2 h 2 
1 + — + — + — + — + — 
8A 8A 2A 2A 6A 



The result in eq. (^) is then recovered, pointing towards restoration in an 
inescapable way. 

Finally, we have checked all the Ward identities in eq. (|T9D ; the explicit results 
for the diagrams involved can be found in Appendix C. 
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Figure 3: One loop self-energy diagrams for the Z gauge boson in the minimal 
standard model. 



3.3.2 Two Higgs doublets 

Models with a richer Higgs structure have several scalar couplings. In order to 
explore symmetry non-restoration, the rule of the game is then to play with the 
freedom in the sign of some of those couplings, while respecting the boundedness 
condition. 

The simplest extension, i.e., the standard model with two Higgs doublets, 
is considered now. We make the usual assumption that the down quarks and 
charged leptons only couple to the Higgs doublet 0i and the up quarks to 02, en- 
suring tree-level flavor conservation of scalar mediated neutral currents. In order 
to avoid radiatively induced FCNC terms, we also impose the discrete symmetry 
0i —0i- The most general, renormalizable, scalar potential consistent with 
the above symmetry and with gauge invariance is: 

V(01, 2 ) = -m^A - ^20202 + Al(0l0l) 2 + A 2 (02> 2 ) 2 

+A 3 (0 t 1 0i)(0 2 > 2 ) + A4I0I02I 2 + ^[A 5 (0 t 1 2 ) 2 + h.c.}. (25) 
The condition for the potential to be bounded from below leads to the constraints: 

X 1 > 0, A 2 > 0, 4AiA 2 > A3, 

4AiA 2 > (A 3 + A 4 + A 5 ) 2 {for A 5 < 0). (26) 
The leading one-loop thermal corrections give the following thermal masses 
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for the fields 4>i,4>2'- 



AVn 



g [ (6Al + 2A 3 + A 4 + ^ 2 + \g' 2 + 3h 2 b + hl)^ 2 
+ (6A 2 + 2A 3 + A 4 + ^ 2 + \g' 2 + 3h 2 t )\<f> 2 \ 2 



(27) 



Although the contributions from both fermions and gauge bosons are positive, 
the scalar couplings A 3 and A 4 may be negative, and therefore it is not possible 
to make any a priori statement about the signs of the mass terms above. What 
can be stated is that the stability conditions in eq.(^) do not allow both mass 
terms in eq. (B7I) to be negative. Since 02 receives a large positive contribution 
from the top Yukawa coupling, it is easier to get a negative thermal mass for 
the field <pi. Then, its vacuum expectation value would remain non-zero at high 
temperature and the SU(2) <g> U(l) symmetry would never be restored. 




Figure 4: Stability bound of the tree level potential (solid line) and parameter 
space leading to symmetry non-restoration (above the dashed-dotted line), for 
A2 = 1 (a) and A2 = 2 (b), with A5 = in both. 



According to eq. fl2"T|), m\(T) < requires 



6A1 + 2A 3 + A 4 + ^g 2 + -g' 2 + 3h 2 b + h\ < . 



(28) 
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Notice that the term |g 2 = 0.99 is already of order one at the electroweak scale, 
which makes it difficult to attain m\(T) < within the perturbative regime 0. 
We have checked numerically that the condition ( pjf ) and the stability bounds 
( |26|) are incompatible for scalar couplings in the range [— 1 < A,; < 1], where 
the weak coupling expression (|27| ) is justified. As an example, in Fig. 4 we plot 
the stability bound of the tree level potential (the allowed range is below the 
solid line) and the curve corresponding to m\(T) = (symmetry non-restoration 
occurs above the dashed-dotted line), for A 5 = and A 2 = 1, 2. 

Heading outside the above mentioned range, the numerical results taken at 
face value seem to point towards the possibility of symmetry non-restoration as 
the scalar sector enters the non-perturbative regime, as can be seen in Fig. 4b. 
Of course, the above computation is meaningless outside that range. 

On the above, we have used as "flag" for symmetry (non-) restoration the 
negative scalar "masses" , that is, the location of the minimum of the potential. 
As discussed in previous section, one could calculate instead the induced gauge 
boson magnetic masses, as an alternative analysis of the fate of the symmetry. In 
the linear realization of the symmetry breaking sector of the minimal standard 
model, the magnetic mass has been computed and proved to show a tendency 
towards vanishing. The temperature corrections to the vacuum expectation value 
of the Higgs field indirectly computed through this procedure agree with the result 
obtained from the effective potential approach. The magnetic mass squared is 
given by g 2 v ^\ In the two doublet case, the magnetic mass squared would be 

given by g 2 " 1 ^ T ^ 2 £2 ; thus if there is a region of parameter space for which one of 
the vev's remains non-zero the magnetic mass will not show a tendency towards 
vanishing and symmetry non-restoration becomes possible. 

From our study of the two doublet model, we conclude that the requirement 
of the validity of perturbation theory points towards the usual assumption of 
restoration of the SU(2) ® U(l) gauge symmetry. In the next section, we extend 
the analysis outside the perturbative regime. 

3.4 Strongly Interacting Higgs Sector 

We now study the behavior of the SU{2) L £g> U{1)y symmetry in the Stan- 
dard Model with a strongly coupled Higgs sector. Strong coupling implies (at 
least naively) heavy physical scalar particles, which can be effectively removed 
from the physical low-energy spectrum. An effective Lagrangian approach is the 
natural technique to use when all the physical degrees of freedom in the sym- 
metry breaking sector are heavy. We then consider the most general effective 
Lagrangian which employs a non-linear realization of the spontaneously broken 
SU{2)l ® U(1)y gauge symmetry p5fl . The resulting chiral Lagrangian is a non- 
renormalizable non-linear a model coupled in a gauge invariant way to the Yang- 

5 We neglect small finite temperature renormalization of the couplings. 
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Mills theory. Chiral Lagrangians have been widely used in the last few years as 
low energy effective theories for electroweak interactions ]34] . 

The Lagrangian keeps only the light degrees of freedom, namely the gauge and 
Goldstone bosons. The latter are collected in a unitary matrix U = exp(iir a T a /v), 
where v is the vacuum expectation value that gives the W and Z gauge bosons 
a mass, n a are the would-be Goldstone fields and r a the Pauli matrices. 

3.4.1 The lowest order Lagrangian 

The lowest order terms in a derivative expansion of the effective Lagrangian are 

v 2 

CochL = -Ti[D^D»U} + C YM + Cgf + Cfp, (29) 



where 



D^U = d^U + i 9 -{W^t)U - i 9 -U{B^) . (30) 



Cym is the pure Yang-Mills piece 



Cym = - l -Tr{W, u Wn - -B^B*", (31) 
and we consider the following gauge-fixing term 

CaF= -\ (w^ wr - J 2 - \ (~k d ^ - 9 '-^ 3 ) 2 ■ (32) 

from which the Faddeev- Popov term, Cfp, can be computed in the usual way. 
The relevant part for our calculation is given in Appendix B. At tree level we 
take £i = £ 2 , so that the gauge boson - Goldstone boson mixing term is canceled. 

Expanding U, we obtain the interaction vertices, in particular the tree level 
masses are given by 

MI = G? 2 + S' 2 )j , M* w = g*f-, (33) 
m^o =£ 2 Mf , ml ± =&M 2 W} (34) 



m2, = V£i6M£ , m c 2 ± = v /£i6M^, (35) 

where n°, ^ are the longitudinal components of the gauge bosons Z, W ± , respec- 
tively, and c°, c ± are the corresponding ghost fields. 

We have computed the one loop temperature corrections to the effective La- 
grangian (|29D at leading order, i.e., 0(T 2 ). Before entering the discussion of the 
results, it is worth to recall the range of validity of the calculation and the approx- 
imations involved. Unitarity implies that this low-energy effective theory should 
be valid for an energy scale much smaller than Attv ~ 3 TeV. Furthermore, as we 
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Figure 5: One loop self-energy diagrams for the gauge bosons in the stan- 
dard model with strongly coupled Higgs sector. Solid lines represent would- 
be-Goldstone bosons, wavy lines gauge bosons and dashed lines Faddeev-Popov 
ghosts. 

shall see in eq.(|36]), our approximation cannot be valid unless T < ^/6v, where 
the latter limit would give the naively extrapolated critical temperature. In the 
vicinity of it, the loop expansion performed here is not appropriate. Therefore, 
our conclusions will be reliable up to T ~ 200 GeV, and expected to be an accept- 
able guideline up to 500 GeV. In order to obtain analytic expressions we work in 
the limit T ^> m,, with rrii the masses of the low energy spectrum (which means 
T ^> gv), and T 3> k, where k are the external momenta. This approximation is 
known in the literature as the hard thermal loop approximation (HTL) [pOt| . 

We are doing an expansion in both T/v and the small coupling constants g, g'. 
The corrections of order g 2 ^ are smaller than which will only appear at 
higher order in perturbation theory since, in the HTL approximation, we assume 
gv <C T. Indeed, notice that already at T = the gauge coupling constants, 
g, g', are not corrected by quadratically divergent diagrams, as can be seen from 
simple power counting arguments on the one-loop diagrams contributing to vertex 
functions: only logarithmic divergences appear. Hence no T 2 renormalization of 
g, g' may appear. 

Let's start with the thermal corrections to the gauge boson masses, obtained 
from the corresponding self-energy tensor (Fig. ^]). It is well known that the 
magnetic mass of the gauge bosons in an unbroken gauge theory vanishes at one 
loop |HJ, so that only diagrams involving would-be-Goldstone boson loops will 
give a non-zero contribution to the magnetic masses in the broken phase, and we 
find 6 

M 2 - q 2V -(i-Ml) 

1VI W,mag " 9 4 I 1 Uy 2 j ' 

6 Although we have not included fermions in our calculation, it is easy to see that they do 
not contribute to the magnetic masses either, at leading order. 
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zj,mag 

Ml 



0. 



where N 



f 



4 \ 12v 2 

2 in SU(2). For the electric masses, defined as ITl(0, k), we get Q 
- 2 / 17 N f T 2 \ 



(36) 



12v 2 



M 



I2\ 



Z,el 



(g 2 + g 



N f T 2 



I2v 2 



NfT 2 
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K9 2 + g' 2 )(c 



.2 
-VI' 



,2 2 



+ og c wi 
(37) 



where cyy {sw) is the cosine (sine) of the weak mixing angle at zero temperature. 
Focusing on the magnetic masses, we can then rewrite them as 



M 2 

lvl W,mag 



9' 



v{Tf 



M 



Z.rnag 



(g 2 + g 



12 



v(Tf 



(38) 
(39) 



with v(T) 2 given by 



v(Tf 



(N- 1)T 2 



12v 2 



(40) 



where N = 3 in SU{2). 

It is worth to remark that the would-be-Goldstone boson field renormalization 
is the same as the one for the Goldstone bosons in the global case (eq. (171)), while 
the temperature corrections to v 2 coincide with those of F£ 2 . 

The diagram in Fig. ^ generates a gauge boson - Goldstone boson mixing 
term proportional to T 2 , which is absorbed by a renormalization of the gauge 
fixing parameter £ 2 ; 

' 2T 2 \ 



6(T) =61 + 



9v 2 



(41) 



With respect to the remaining parameter in the effective Lagrangian in eq. 
), the gauge fixing one £i, it is not renormalized at 0(T 2 ). The same applies to 
the gauge boson and ghost fields. Naive dimensional counting shows that the 1PI 
one-loop diagrams that will renormalize those entities are at most logarithmically 
divergent, alike to the situation for g and g' discussed earlier, and thus not able 
to produce T 2 corrections. 



We acknowledge C. Manuel for pointing out two misprints in these formulae, in an earlier 
version. 
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Figure 6: One loop diagram which generates the 7r — W mixing term. 




Figure 7: One loop self-energy diagram for the Faddeev-Popov ghosts. 



The consistency of our results has been verified by computing the leading 
order corrections to the masses of ghost and would-be-Goldstone bosons through 
the corresponding one loop self-energies (Figs. [7| and ||). Once the would-be- 
Goldstone boson field renormalization has been taken into account, the temperature- 
dependent masses depend on the renormalized parameters in the same way as at 
tree level, i.e., 

ml ± {T) = & {T)M*r t mag , (42) 
m 2 AT) = ^6(T)M^, (43) 

The one-loop effective Lagrangian does not have exactly the same functional 
form as the original bare one, eq. (29): it splits, as exemplified by the differing 
electric and magnetic masses. Eqs. (40) and (41) allow to connect several impor- 
tant finite T quantities in a compact notation, though. Eqs. (42) and (43) are 
an example of it. 

We have also checked all the Ward identities in eq. ([H]). The explicit results 
for the different diagrams involved are given in Appendix C. 

A natural question is the relationship with the linear case discussed in subsect. 



3.3. 1[ one expects that taking there the Higgs mass to infinity the results of the 
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Figure 8: One loop self-energy diagrams for the would-be-Goldstone bosons. 



present section should be recovered. A superficial look does not show it. For 
instance, taking the limit A — > oo in eq. fl25|) for the magnetic mass in the 
linear case, eq. (|3q) is not recovered. There is no inconsistency, though: the 
high temperature and heavy Higgs mass limits are not interchangeable. We have 
indeed checked that eq. (j36|) is obtained by taking the limit m — > oo (m being 
the physical Higgs mass), before doing the one-loop computation in the linear 
case (which implies not considering the following diagrams of Fig. 3: 2, 5, 10, 11 
and 7, 9 when the physical Higgs boson is in the loop). The same argumentation 
is valid for the rest of the physical parameters. 

Regarding the question of the SU(2)<S>U(1) symmetry non- restoration which 
initially motivated our study, we conclude from eq. fl36|) that thermal effects 
tend to restore the symmetry also in the non-perturbative regime. Notice that 
although the magnetic mass is non-vanishing in the symmetric phase beyond one- 
loop order, it is expected to be of order g 2 T, and therefore much smaller than 
the magnetic mass in the broken phase, of order gv (recall that g 2 T <C gT < gv). 
Thus we can interpret the decreasing of the magnetic mass with the temperature 
as a "flag" for symmetry restoration. On the contrary, for instance the W electric 
mass in eq. fl37|) can be written as 

n,2 2V 2 (T) 2 3N f T 2 

M^ el = g 2 -f 1 + g 2 — |— , (44) 

and in the symmetric phase, while v(T) = 0, it is nevertheless non-vanishing and 
of order gT, as anticipated. 

These results are shown qualitatively in Fig. 9. We plot both the magnetic 
and the electric masses of the W gauge boson as functions of the temperature, 
at leading order T 2 . The solid and dashed lines correspond to our one loop 
calculation, and the dotted line to the non-perturbative estimate of the magnetic 
mass in the symmetric phase, M^™ a9 = 0.28g 2 T, which is taken from ref. [IS . 
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Figure 9: Electric mass (dashed line) and magnetic mass (solid line) for the W 
gauge boson in the broken phase, and non-perturbative estimate of the magnetic 
mass (dotted line) in the symmetric phase, M^-™ afl = 0.28g 2 T. 



3.4.2 Model dependence 

As already mentioned, the lowest order term in the derivative expansion of the 
effective Lagrangian, CcchL, has a universal character. The next term in the 
expansion, is model dependent, namely, it depends on the specific dynam- 
ics of the symmetry breaking sector through the different values of the various 
constants. The logarithmic divergences generated at one loop by Cgchl are con- 
sistently absorbed by the renormalization of those constants. As stated before, 
we have neglected all the zero temperature renormalization effects, and therefore 
the model dependence contained on them, as well as in the matching conditions 
0. This is justified, since we are looking for temperature effects. 

The only model dependent contribution to the pure thermal corrections at 
one loop is due to the dimension two term 

C p = \pv\Tr{Ur z U\D,U)U^))\ (45) 

which explicitly breaks the custodial SU(2)c symmetry. This term contributes 
to Ap at tree level, and is thus strongly constrained by experimental data. The 

8 We thank J. Matias for pointing out this fact to us. 
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Figure 10: Contribution of the model dependent term Cp to the one loop self- 
energy of the W and Z gauge bosons. 



contribution of Cp to the magnetic mass of the gauge bosons is given by 

SM^ mag = ~9 2 (3^ (46) 

SMl mag = g 2 ^. (47) 

The sign of the parameter (3 can be either positive or negative. Whatever the 
case for a given theory, eqs. (46) and (47) show an opposite behaviour for the 
W and Z magnetic masses, which are no more forced to behave alike since the 
operator under study breaks the custodial SU(2)c symmetry. When combined 
with the universal leading contribution found in eqs. (|56|), the total correction 
reads 



2 T 2 T 2 



Z ' ma 9 = ^-^ + ^^+ 9 -)v^ ^ 

The SU(2) <S>U(1) symmetry can be considered effectively restored only when 
all possible "flags" have signaled it. The above result, taken at face value, would 
indicate that the SU(2) £g> U(l) symmetry may never be restored at high temper- 
ature for theories where a large enough value of the coefficient of the operator 
C/3 is generated. Such a strong statement has to be tempered by recalling that, 
if the chiral expansion is valid, we expect (3 to be small (in typical models it is 
of order of a coupling constant squared), and the total correction in ( f4Tf ) would 
be dominated by the leading one, pointing in a natural way towards restoration. 
Moreover, low energy constraints on new physics give an experimentally allowed 



value of j3 of order 10 [32], which implies that in phenomenologically acceptable 
models the contribution of the operator Cp is indeed negligible. It is interesting 
to retain that the tendency to symmetry restoration can be reversed for values 
of j3 which are not outrageously large, as seen from eqs. (§^). 
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4 Conclusions 



We show that the spontaneously broken SU{2) <8> £7(1) gauge theory in mod- 
els where the Higgs sector becomes strongly interacting (such as composite Higgs 
models and technicolor-like ones) tends to be restored when the system is heated. 
This conclusion is obtained in a model-independent way using the techniques of 
the electroweak chiral Lagrangian. Specific models will only affect the sharpness 
of such a tendency, unless the natural chiral expansion is not respected. We 
quantify such model dependence computing the generic contribution of the lead- 
ing effective operator whose coefficient is model-sensitive, C@; its one-loop con- 
tribution to the W and Z magnetic masses is found to have opposite sign. The 
technique, while only valid for temperatures lower than the electroweak scale, has 
the advantage of its non-perturbative character. The physical conclusion reached 
here parallels the corresponding one for the other main avenue of beyond the 
standard model physics, supersymmetry, where perturbative treatments show a 
tendency towards restoration. 

In this work we have as well explored the £77(2) <S>U(1) gauge symmetry in 
a perturbative regime: the cases of one and two light Higgs doublets. Again, 
the results show symmetry restoration at high temperatures when the full scalar, 
gauge boson and fermion corrections are taken into account. 

The above conclusions have been obtained mainly through the study of the 
temperature dependent magnetic mass for the gauge bosons, which we propose 
as an appropriate "flag" in the broken phase. 

In addition, BRST invariance has been explicitly checked for gauge theories 
at finite temperature, a novel result. 

Finally, it is worth to remark that global symmetries have been studied as 
well for the non-linear sigma model at finite temperature. While this subject 
and the results are not new, the technical approach we used is so: we derive first 
the temperature corrected one-loop effective Lagrangian, from which the physical 
conclusions are then extracted. 

A grain of salt: we have disregarded the putative role of finite parameters 
such as finite volumes or causal domain sizes in the history of the universe. Their 
effect could constitute an interesting topic to study. 

Appendix A: Tensor Basis 

The tensor basis in terms of which we have expressed the gauge boson self-energy 
is given by, 



(50) 




(51) 
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K fl K u + K^K V 
&v = a a + a a (52) 

= — (53) 

where A' 1 = [K ■ uK^ - K 2 vf) /k and fc is such that A^A^ = cj 2 - A; 2 with 
w = A^u". 

Appendix B: Faddeev- Popov Lagrangian 

The Faddeev- Popov Lagrangian which corresponds to the non-linear realization of 
the S77(2)(g){7(l) gauge symmetry is different from the one derived for the minimal 
standard model for which the gauge symmetry is linearly realized. Here we 
present the Faddeev- Popov Lagrangian terms which are relevant for our purposes. 



More general results can be found in [33 . 



C FP = c+[-V 2 - (^)g - + tt 2 ) + ...]]c + (54) 

2 2 bv 

E ct[-V 2 -( 9 ^)[l-l^ + nl) + ...]]c l + 

(cfc 2 -4cd[-gd»W* + ( 9 -^)^} + 
(ctcs-cic l )[gd^ + ( g -^)^] + 
(cic 3 -4c 2 )[gd^ + ( 9 -^)^] + 
gg'^j(clc 3 + c\c a )[v - ^-(tt 2 + tt 2 ) + ...] + ... 



where £i = £ 2 = f • 

Appendix C: Ward Identities at One Loop 

We have verified the Ward identities at one loop, both for the linear and non-linear 
realization of the gauge symmetry, for the two-point Green functions, computing 
the W^, Zp, gauge boson, Goldstone boson and ghost self-energies, together 
with the one loop Goldstone boson-gauge boson mixing term. At leading order 
(0(T 2 )) and for small external momenta the results for the standard model case 
are: 



d vy ) g Ky ^ y J Sm 2 ' 
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T 2 1M 2 T 2 

yrW _ 2 1 n 2 01VI W 1 

iiD " 9 T~ 9 ^m 2 ~> 
W D = 0, 

up = 0, 

IT = 0, 



16M W * 16m 2 ' 



where m 2 = 2/i 2 represents the Higgs boson mass squared. 

Concerning the electroweak chiral Lagrangian the results for the two point 
Green functions are: 





= i /r • g 


n D 


rp2 
2 1 

1 ^24' 


TT 7 


= o, 


n D 


= o, 


W 


rp2 

~ 18t; 2 ' 


IT 7 


= o, 




T 2 

= t9 lBv> 



rp2 

/2\_ 

J 24' 



n c± = -ig 2 ^. (56) 
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